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^ ■ Abstract 

o ■ We study a Dirac neutrino mass model of Davidson and Logan. In the model, the smallness of 

,^ , the neutrino mass is originated from the small vacuum expectation value of the second Higgs of 

^- _ two Higgs doublets. We study the one loop effective potential of the Higgs sector and examine how 

> 

\^ • the small vacuum expectation is stable under the radiative correction. By deriving formulae of the 

f^ . radiative correction, we numerically study how large the one loop correction is and show how it 



depends on the quadratic mass terms and quartic couplings of the Higgs potential. The correction 
changes depending on the various scenarios for extra Higgs mass spectrum. 

PACS numbers: 12.60.Fr,14.60.St,14.80.Ec,14.80.Fd 



I. INTRODUCTION 



The smallness of the neutrino mass compared with the other quarks and leptons are one 
of the mysteries of nature. Recently, a new mechanism generating small Dirac mass terms 
for neutrino has been proposed l|-l3|. The similar mechanism generating the small neutrino 
Dirac mass term for the TeV seesaw mechanism is also proposed in [-^ and phenomenology is 
studied in |5| and |6|. There are also models with radiatively generated Dirac mass term in 
71, l8|. The interesting feature of the modelproposed in [l|, |2| is the tiny vacuum expectation 
value for an extra Higgs SU(2) doublet [9|. The small neutrino mass is realized without 
introducing tiny Yukawa coupling for neutrinos. A softly broken global U(l) symmetry 
guarantees the tiny vacuum expectation value for the extra doublet. In addition to the 
small softly breaking mass parameter, the mass squared parameter for the extra Higgs is 
chosen to be positive so that the light pseudo Nambu Goldstone bosons due to the softly 
broken global symmetry do not appear. This is a contrast to the mass squared parameter 
for the standard model like Higgs boson. 

In the present paper, we study the global minimum of the tree level Higgs potential by 
explicitly solving the stationary conditions. There are many studies of the tree level Higgs 
potential of general two Higgs doublet model 10l-ll5|. (See also 16[ for recent review of 



two Higgs doublet model.) It has been shown that the charge neutral vacuum is lower than 
the charge breaking vacuum [lO| . Also the vacuum energy difference of two neutral minima 



was derived 



12 



14| . We make use of the results and identify the vacuum of the present 



model. When the U(l) symmetry breaking term is turned off, the tree level Higgs potential 
and the phase structure of the present model is rather similar to the model with Z2 discrete 
symmetry [l7|,|l8[. In contrast to Z2 symmetric case, it is essential to keep the soft breaking 
term when finding the true vacuum. If we set the symmetry breaking term at zero, then 
the order parameter corresponding to the softly broken U(l) symmetry becomes redundant 
parameter and can not be determined. We treat the soft breaking term as small expansion 
parameter and obtain the vacuum expectation values and the vacuum energies in terms of 
the parameters of the Higgs potential. 

The constraints on the parameters of the model for which the desired vacuum can be 
realized, are derived and they are rewritten in terms of Higgs masses and a few coupling 
constants which can not be directly related to the Higgs masses. These constraints are fully 



used when we study the radiative corrections to the vacuum expectation values numerically. 

Beyond the tree level, we study the radiative correction to the Higgs potential and the 
vacuum expectation values of Higgs. Since the neutrino masses are proportional to the 
vacuum expectation value of one of Higgs , one can also compute the radiative corrections 
to neutrino masses. As already noted in [l| the radiative correction to the softly breaking 
mass parameter is logarithmically divergent and it is renormalized multiplicatively. We 
derive the formulae for the one loop corrected vacuum expectation values for two Higgs 
doublets by studying one loop corrected effective potential. The corrections are evaluated 
numerically by exploring the parameter regions allowed from the global minimum condition 
for the vacuum. We show how the radiative corrections change depending on the extra Higgs 
spectrum. The radiative corrections are also evaluated for the case that a relation among 
the coupling constants is satisfied. 

The paper is organized as follows. In section (II), we derive the condition for the desired 
vacuum being global minimum. In section (HI) one loop effective potential is derived and 
one loop corrections to the vacuum expectation values are obtained in section (IV). In section 
(V), the corrections are evaluated numerically for various choices of parameters of the Higgs 
potential. Section (VI) is devoted to summary and discussion. 

II. MODEL FOR DIRAC NEUTRINO WITH A TINY VACUUM EXPECTATION 
VALUE 



fl, 



The model of the Dirac neutrino is proposed in [l|. In [ij, two Higgs SU(2) doublets are 
introduced, 

where $i's vacuum expectation value is nearly equal to the electroweak breaking scale and 
the second Higgs $2 has a small vacuum expectation value which gives rise to neutrino mass. 
The Higgs potential in []J is, 

Kree = J] (ml^l^, + ^($1$.)') - {ml,^\^2 + k.C.) + A3($I$l) ($^$2) + XM^2\'.{2) 
i=l,2 ^ ^ 

U(l)' charge is assigned to the second Higgs. The U(l)' global symmetry is broken softly 
with the term m^2- ^^ this paper, we introduce the following real 0(4) representation for each 



doublet, because this parametrization is convenient when computing the one loop corrected 
effective potential. 
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Using the notation above, the tree level effective potential introduced in Eq.([2]) can be 
written as, 
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where one can choose mf^ real and positive. With the notation of Eq.(|3]), the softly broken 
global symmetry U(l)' corresponds to the following transformation on (p2, 
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01 does not transform under U(l)'. Therefore U(l)' is broken softly when m^2 does not 
vanish. Without loss of generality, one can choose the vacuum expectation values of Higgs 
with the form given as. 
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' f sin /3 sin a COS ^' ^ 

—V sin /3 sin a sin 6' 

V sin (3 cos a cos 6' 

\—v sin /3 cos a sin 6' J 



(6) 



where the range for 9' is [0, 2ir) and the range for /3 and a is [0, ^]. We call the four order 
parameters as (fi = {v,/3,a,9'), {I = 1,2,3,4). When mi2 vanishes, by taking = ^' in 



Eq.®, one can rotate 9' away in Eq.([6]). For the most general case, in total, there are four 
independent order parameters when U(l)' symmetry is broken. 

For completeness of our discussion, we give the constraints on the quartic couplings from 
condition that the tree level potential is the bounded below, [l|, 10|, 19|. 



Ai > 0, A2 > 0, (7) 

-\A^<A3, (8) 

-\/V^<A3 + A4. (9) 

In addition to the conditions on the quartic terms, one can constrain the parameters includ- 
ing the quadratic terms so that the desired vacuum satisfies the global minimum conditions 
of the potential. About the global minimum of the tree potential, it was shown that the 
energy of charge neutral vacuum is lower than that of the charge breaking vacuum 10| . We 
therefore set a zero. We also require the vacuum expectation value of the second Higgs is 
much smaller than that of the first Higgs, which implies that tan j3 is small. In terms of the 
parametrization in Eq.(|6]) with a = 0, the potential can be written as, 

V,,,,{v, /?, 6') = A{l3y + 5(/3, e')v\ (10) 

where, 

Al „_ 4 n I A2 . 4 „ / A3 A4\ 2 a ■ 2 



A{P) = -^ cos^ /3 + ^ sin^ /3 + I _i + ^ 1 cos^ /? sin^ /3, 

2 2 

B{(3,9') = -^ cos^ f3 + ^ sin^ 13 -ml^ cos 6' cos 13 sin /3. (11) 

We first find the global minimum of Kree- The stationary conditions ^ '""'' = (/ = 1, 2, 4), 
are written as, 



-') 



v{2Av^ + B) = 0, (12) 

o • no (l-rir2)cos2^ + r2-rir3 

2r4 = sm2n — ; ^ , (Id) 

^r2Cos2 2/3 + (r3 + l)cos2/3 + r2 ^ ^ 

m^2sin6''sin2/3 = 0, (14) 



where rJi 



4) are defined as, 
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Ai + A2 - 2A3 
ml2 cos 9' 

"^11 + "^22 ' 
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The stationary conditions Eq. (TT2|) and Eq. (fT3|) correspond to Eq.(36) of [141]. Here we 
solve them explicitly by treating the soft breaking term rriu as perturbation. The non-zero 
solution for f^ in Eq. (fT2|) is written as, 
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mli + 77122 1 + '"i COS 2/3 — 2r4 sin 2/3 



where A 
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Ai + A2 - 2A34 cos2 2/3 + r3 + 2r2 cos 2f3 
A3 + A4. Substituting it into Kree, one obtains. 
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2(Ai + A2 - 2A34) cos2 2^ + 2r2 cos 2/3 + r3 
For non-zero 771^2 and sin 2/3, the solution of Eq.l TT^ is sin^' = 0. One still needs to find /3 
among the solutions of Eq. (IT3|) . which leads to the minimum of Knin.- We solve Eq.f fT3|) and 
determine /3 by treating r^ (m^2) ^s a small expansion parameter. One can easily find the 
approximate solutions as. 
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cos 6*' = sign(m^^A2 - m|2A34), 
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(3) COS2/3 = ^L(A34+A.)-mi,(A34+A,) ^^^ 

V. ^ '' "^ mf^(-A34+A2)+m.^2(--*'34+-^l) ^ ^ 

Corresponding to each solution, (1)~(3) of Eq.f fTSj) . the vacuum expectation value v"^ and 
the minimum of the potential are obtained. 
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TABLE I: Classification of tlie solutions witli non zero sin 2/3 of the stationary conditions of Higgs 
potential. For (3), 0(r4) correction is not shown. 
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TABLE IL Classification of the solutions with sin 2/3 = 0. 

The leading terms of the vacuum expectation values agree with those obtained in Z2 sym- 
metric model 18|. If sin 2/3 = 0, then r^ must be vanishing and cos^' = from Eq. (fT3|l and 



Eq. flMl) . The vacuum energies of the non-zero sin 2/3 solutions are shown in Tables [H In 
Table [TTl the vacuum energies of the solutions with sin 2/3 = are summarized. 

Next we derive the constraints on the parameters so that the solution corresponding to 
(1) in Table [l] becomes the global minimum of the potential. Since the other cases (2)-(5) 
do not have desired properties , we restrict the parameter space so that these solutions can 
not be a global minimum. Since v must have large positive vacuum expectation value, m^^^ 
must be negative. In order that the vacuum energy of (1) is lower than that of (4) , 

ml^Xi - ml^Xu > 0, (cos^' = 1). (20) 

When Eq. (!20!) is satisfied and the solution (1) does exist, one can show that the vacuum 
energy of solution (3) is higher than that of (1). Furthermore when 17122 > O5 the solutions 
corresponding to (2) and (5) are not realized. Then one can state the region of parameter 
space which is consistent with the case that the vacuum (1) becomes global minimum is. 
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m 



m 



22 \ 
11 



(21) 



Next we consider the case with negative 77122 • I^ this case we impose the additional condition 
so that the vacuum energies corresponding to (2) and (5) are higher than that of (1). 
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Then the condition for (1) is global minimum in this case is, 
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m^i < 0, m^2 < 0, A34 > — pAi, Aa^ > Ai— |^. (23) 



In the following sections, we explore the regions for the parameters obtained in Eq. lpTI) . 
Eq.(l23]), Eq.® and Eq.©. 

III. EFFECTIVE POTENTIAL IN ONE LOOP AND RENORMALIZATION 

In this section, we derive the effective potential within one loop approximation. We in- 
troduce a real scalar fields with eight components as, 0* = (0j, (jyf, (pl, (pf, (pl, (/)l, (pl, (j)2Y ,{i = 
1 ~ 8). With the notation above, the one loop effective action is given as, 

rT" = ^llndetD-\cP), 
D-^ = n + M|, (24) 

where M^ is the mass squared matrix of the Higgs potential, 

„ „ I mi, X 1 \ „ 

M^ = M\(l>)+\ \-ml,ai, 

\ 171^2 y< ly 

d(f)id(t)j 
where 1(0) denotes 4x4 unit (zero) matrix, ai is defined as. 

In Eq.( l26l) . 1(0) also denotes a four by four unit (zero) matrix. In modified minimal subtrac- 
tion scheme, the finite part of the one loop effective potential becomes, 

A-d r ^d^ 
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' Tr. f M^(LnM _ 3) , (27) 
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Vc denotes the counterterms and the derivation of Vc can be found in Appendix A. 



IV. ONE LOOP CORRECTIONS TO THE VACUUM EXPECTATION VALUES 



In this section, we compute the one loop corrections to the vacuum expectation values. 
Using the symmetry of the model, in general, one can choose Lpi = [v, (3, a, 6') as the vacuum 
expectation values of Higgs potential. Their values are obtained as the stationary points of 
the one loop corrected effective potential V = Vtree + Viioop, 

dV 



dipi 



0. 



(28) 



By denoting the vacuum expectation values as sum of the tree level ones and the one loop 



corrections to them; 9?/ 
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where Mr, is a diagonal 8x8 tree level mass squared matrix of Higgs sector and Ljj is 
4x4 matrix given by the second derivatives of the tree level Higgs potential with respect 
to the order parameters. 
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The diagonal Higgs mass matrix squared M|, is related to 8 x 8 Higgs mass matrix squared 
M| in Eq.(l25]). 
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where M|-q is obtained by substituting the vacuum expectation values to M^. O is shown in 
appendix D. Since Mr, is the 8x8 diagonal matrix which elements correspond to the Higgs 



masses and zero mass of the would be Nambu-Goldstone bosons, one may write Eg. (12 9 p in 
a simple form. The Higgs masses squared in Eq. fl3T]) are given by, 
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(32) 



where 7 is an angle with which one can diagonalize the 2x2 mass matrix for CP even 
neutral Higgs. tan 27 is given as, 

-4m22 + 2sin2/?(A3 + A4)i'2 



tan 27 



(33) 



(3(-Ai cos2 /3 + A2 sin^ /3) + cos 2/3(A3 + A4))f ^ - 2(ml^ - ml^) ' 

To compute Eq. fl29|) . we still need to calculate O^^^O and Ljj. They are shown in ap- 
pendix C. Using Eq.( l29|) and Eq.( ICl|) . one can find the quantum corrections for a and 6' 
vanish, 
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For v^^^ and I3^^\ one obtains. 
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where Ll is 



i'=r"^'M. (36) 

The elements of L' are shown in Eq. flC4l) . Eg. (1351) corresponds to the one loop exact formulae 
and is a main result of present paper. In the leading order of the expansion with respect to 
the symmetry breaking term mfj, the correction to v becomes, 

The Higgs masses in the formulae are the ones in the limit of m\2 -^ 0, 

Ml ~ m\^ + ^Ai^;2, 
Ml ^ Ml c. ml, + h±h,^, 

Ml, ^ ml, + ^v\ (38) 





where v is related to mf^^ as, 

^V^-m?!. (39) 

The approximate formulae for the physical Higgs masses in Eq. ( l38l) which are valid the limit 
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2 — ;■ 0, agree with the ones given in p| except the notational difference of Mh and M^ 



2l| . The one loop correction to (5 in the leading order expansion of m^g is given as. 
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(41) 



Eq.( HO|) shows that the quantum correction is also proportional to the soft breaking param- 
eter 771^2 which is expected. We also note that the correction depends on the Higgs mass 
spectrum and quartic couplings. The correlation to Higgs spectrum is studied in the next 
section. 

V. NUMERICAL CALCULATION 

In this section, we study the quantum correction to j3 and v numerically. As shown in 
Eq. fl371) and Eq. fHOl) . the quantum corrections are written with four Higgs masses and the 
four quartic couplings. Since the neutral CP even and CP odd Higgs of the second Higgs 
doublet are degenerate as Ma = Mh in the limit rriu -> (See Eq.(l38l)). the three Higgs 
masses {Mh,Ma,Mh+) are independent. Moreover for a given charged Higgs mass and 
neutral Higgs mass, Ai and A4 are given as, 

'^l — — 5-1 

Ml - Ml. , , 

A4 = 2 ^ , ^+ . 42 

yA 

A2 and A3 are the remaining parameters to be fixed. The lower limit of A3 obtained from 
Eq.dH]) and Eq.(l9]) is written as, 

Ma..(-^V/A;,-^V^-2MZ^,<A, (43) 

One can also write A3 with the charged Higgs mass formulae, 

\^ = \{Ml^-m\^). (44) 

Depending on the sign of ni^^, the upper bound and the lower bound of A3 can be obtained 
for a given charged Higgs mass. Combining it with Eq.( H3l) . the constraints for positive 777-22 
case are. 

Max. ^v^, -v^-2 -^ , ^^ <A3< 1^, m^2>0- 45 

When 77^22 < 0, in addition to the lower bound on A3, the constraint on A2 in Eq. (l22!) should 
be satisfied, 

^ < A,, v/a; XA3 - 2^)-^. «<0)^ (46) 
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0.14 
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0.42 
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0.55 


0.69 
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TABLE III: The coupling constants (A3,A2) which satisfy the relation, Eq. (|47|) for the three de- 
generate masses Mh+ = Ma = 100,200 and 500 (GeV). 

Now we study the quantum corrections numerically. We fix the standard model like Higgs 
mass as Mh = 130 (GeV). There are still four parameters to be fixed and they are A2, A3, Ma 
and Mh+ ■ Focusing on the Higgs mass spectrum of the extra Higgs, we study the radiative 
corrections for the following scenarios for Higgs spectrum and the coupling constants. 

A. Case for Ma = Mu+; degenerate charged Higgs and pseudoscalar Higgs and a 
relation for vanishing quantum correction f3^^' 

We first study the corrections for degenerate charged Higgs and pseudoscalar Higgs. In 
this case, for a given degenerate mass, one can identify the values of coupling constants A2 
and A3 for which (3^^'^ vanish. With Ma = Mh+^ the relation for coupling constants which 
satisfies /3^^^ = is, 
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(47) 



The set of coupling constants (A3, A2) which satisfy the relation Eq.f H7|) are shown in table 
HI. We note that when A2 is as large as 10, A3 is at most about 3. If A2 is 1, A3 is lies in the 
range 0.55 ~ 0.7. 
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FIG. 1: The quantum correction ^-g- (gray lines) and ^^ (black lines) due to the non-degeneracy of 
charged Higgs and pseudoscalar Higgs masses. The pseudoscalar Higgs mass Ma (GeV) dependence 
of the quantum corrections - — {x = /3, v) is shown while the charged Higgs mass is fixed as 
Mff+ = 100 (GeV). The set of parameters (A3, A2) are chosen so that the correction I3^^> vanishes 
for the degenerate case; M^+ = Ma = 100 (GeV). The values (A3,A2) are taken from Table (III) 
and they are (0.19,0.14) (solid line), (0.28,0.28) (dashed line), (0.41,0.56) (dotted line), (0.55,1) 
(dotdashed line), and (1.8,10) (thick solid line). 

B. Non-Degenerate case Ma / Mu+ with the coupling constants satisfying Eq. ()47p 

Next we lift the degeneracy by shifting the pseudoscalar Higgs mass from the charged 
Higgs mass and study the effect on [5^^^ and v^^\ The non- degeneracy of the charged Higgs 
mass and the pseudoscalar Higgs mass is constrained by p parameter. We change the 
pseudoscalar Higgs mass within the range \Ma — Mh+\ < 100 (GeV) allowed from the 
electro-weak precision studies. The coupling constants (A3,A2) are chosen from the sets of 
their values satisfying the relation Eq. fH7|) . In Fig.l, we show ^^ as a function of Ma with 
charged Higgs mass Mh+ = 100 (GeV). When Ma = 100 (GeV), the correction vanishes 
exactly. As we increase Ma from 100 (GeV) (the mass of charged Higgs), the correction 
becomes non-zero and is negative. The corrections are at most about 1.3% when A2 ~ 1. 
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FIG. 2: The quantum correction ^-g- (gray lines) and ^^ (black lines) due to the non-degeneracy of 
charged Higgs and pseudoscalar Higgs masses. The pseudoscalar Higgs mass Ma (GeV) dependence 
of the quantum corrections - — (x = /3, v) is shown while charged Higgs mass is fixed as Mij+ = 200 
(GeV). The set of parameters (A3,A2) are chosen so that the correction ji^"^' vanishes for the 
degenerate case; Mu+ = Ma = 200 (GeV). The values (A3, A2) are taken from Table (III) and they 
are (0.16,0.14) (solid line), (0.28,0.28) (dashed line) , (0.47,0.56) (dotted line), (0.69,1) (dotdashed 
line), and (2.8,10) (thick solid line). 

By increasing Ma further, we meet the point around at Ma — 200 (GeV) corresponding 
to that the correction vanishes again. In Fig. 2, we study the correction (3^-^^^ with larger 
charged Higgs mass case, Mh+ = 200(GeV). In contrast to the case for Mh+ = 100 (GeV), 
by increasing Ma from 200 (GeV) where the correction vanishes, it increases and becomes 
positive. We also note that the correction tend to be larger than the lighter charged Higgs 
mass case. When A2 ~ 1, increasing the pseudoscalar Higgs mass from 200 (GeV) to 300 
(GeV), the correction is about 10%. As the pseudoscalar Higgs mass decreases from 200 
(GeV) to 100 (GeV), the correction becomes negative for < A2 < 1. With the larger value 
A2 = 10, we meet the point around at Ma — 150(GeV) where the correction vanishes again. 
In Fig. 3, we study the further larger charged Higgs mass case, i.e., Mh+ = 500 (GeV). With 
Ma — 600 (GeV), the correction is positive and about 100%. The correction stays small for 
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x("/x 




Ma (GeV) 



FIG. 3: The quantum correction ^o- due to the non-degeneracy of charged Higgs and pseudoscalar 
Higgs masses. The pseudoscalar Higgs mass M/i(GeV) dependence of the quantum corrections - — 
(x = /3, v) is shown while charged Higgs mass is fixed as Mu+ = 500 (GeV). The set of parameters 
(-^3)'^2) Eire chosen so that the correction P^^' vanishes for the degenerate case; Mij+ = Mj^ = 
500 (GeV). The values (A3,A2) are taken from Table (HI) and they are (0.18, 0.14) (solid line), 
(0.28,0.28) (dashed line), (0.42,0.56) (dotted line), (0.59,1) (dotdashed line), and (2, 10)(thick 
solid line). 



< A2 < 1 when decreasing Ma from 500 (GeV) to 400 (GeV). 



C. The correction 



„{i) 



In Figures 1,2, and 3, we also show the correction - — as functions of Ma- v^^"* is 
independent on A2 and does not necessarily vanish at the same points where /J'-^^ vanishes. 
With A3 > 2 and Mh+ > 200 (GeV), when the pseudoscalar Higgs mass is much larger 
than that of charged Higgs mass, we find very large correction to v. In Fig. 4, we show that 
the two dimensional surface which corresponds to v*^^-* = 0. We find that the interior of 
the surface corresponds to the region of the positive correction; v^^^ > while the exterior 
region of the surface corresponds to the negative correction; v^^' < 0. 
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FIG. 4: The two dimensional surface for v^^' = 0. 




FIG. 5: The regions of {Mjj+,Ma) which cor- 



respond to 



,3(1) 



(0,0) (dark gray), 



(0.01,0.01) (gray), and (0.1,0.1) (hght gray). 



In Fig.5, we have shown the regions of {Mfj+,MA) which correspond to that the cor- 
rections of 1^;^^^! and |/3*-^^| have the definite values (0,0.01,0.1). The dark gray shaded 
area corresponds to the region where both v^^'^ and /3*-^^ can vanish with taking account of 
the conditions; Eq.(I7]), EqdHD and Eq.(l9]). We note that for Mh+,Ma > 200 (GeV), the 
quantum corrections vanish around the region where the charged Higgs degenerates with 
the pseudoscalar Higgs. When the corrections become larger, the larger mass splitting of 
the pseudoscalar Higgs and charged Higgs is allowed. However as the average mass of the 
charged Higgs and pseudoscalar Higgs increases, the allowed mass splitting becomes smaller. 

VI. DISCUSSION AND CONCLUSION 

In this paper, the Dirac neutrino mass model of Davidson and Logan is studied. In the 
model, one of the vacuum expectation values of two Higgs doublets is very small and it 
becomes the origin of the mass of neutrinos. The ratio of the small vacuum expectation 
value V2 and that of the standard like Higgs Vi is tan/3 = — . Therefore tan/3 is very small 
and typically it is 0(10"^). The smallness of tan/3 is guaranteed by the smallness of the 
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soft breaking term of U{iy. 

We have treated the soft breaking term as perturbation and calculated, in particular, 
the vacuum expectation of Higgs in the leading order of the perturbation precisely. As 
summarized in Table I, only by including the soft breaking terms, one can argue which 
of the local minima minimizes the potential and becomes the global minimum. We have 
studied the global minimum of the tree level Higgs potential including the effect of the soft 
breaking term as perturbation. 

Beyond the tree level, we study the quantum correction to the vacuum expectation values 
and tan/3 in a quantitative way. In one loop level, we confirmed that tree level vacuum 
is stable, i.e., the order parameters which vanish at tree level do not have the vacuum 
expectation value as quantum correction. In one loop level, we derived the exact formulae for 
the quantum correction to /3 in the leading order of expansion of the soft breaking parameter 
m^2- We have confirmed not only that the loop correction to tan/3 is proportional to the 
soft breaking term but also found that the correction depends on the Higgs mass spectrum 
and some combination of the quartic coupling constants of the Higgs potential. Technically, 
we carried out the calculation of the one loop effective potential by employing 0(4) real 
representation for SU(2) Higgs doublets. 

Dependence of the corrections on the Higgs spectrum is studied numerically. We first 
derive a relation of the coupling constants which corresponds to the condition that the 
correction to /3 vanishes for degenerate extra Higgs masses. Next, we study the effect 
of non-degeneracy of the charged Higgs and pseudoscalar Higgs on the correction. If the 
charged Higgs mass is as light as 100 (GeV) ~ 200 (GeV), allowing the mass difference of 
charged Higgs and pseudoscalar Higgs is about lOO(GeV), the quantum corrections to both 
(3 and v are within a few % for (A3, A2) ~ (0.5, 1). If the charged Higgs is heavy Mh+ = 500 
(GeV), a slight increase of the pseudoscalar Higgs mass from the degenerate point leads to 
very large corrections to /3 and v. 

One can argue the size of the quantum corrections to the neutrino mass of the model, 
because the ratio of the tree level neutrino mass and one loop correction can be written as, 

^ = — + ^, (48) 

where we take account of the corrections only due to Higgs vacuum expectation values. The 
formulae Eq. fHHl) implies that radiative correction to neutrino mass is related to the Higgs 
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mass spectrum. Therefore once Higgs mass spectrum is measured in LHC, one can compute 
the radiative correction to the mass of neutrinos using the formulae Eq. fHHj) . 
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Note added 

After submitting the pa per, we are aware that the stability of the model studied in this 



paper was also discussed in [20(|. Compared to their analysis, we derived the 1 loop effective 
potential taking into account all the interactions of Higgs sector while they consider a part 
of the interactions and study the stability in a qualitative way. Using the effective potential, 
we carried out the quantitative analysis of the quantum corrections. 

Appendix A: Derivation of one-loop effective potential 

In this appendix, we give the details of the derivation of the one-loop effective potential 
and the counterterm in Eq. fl27l) . One can split M^(0)jj in Eq.f l25p into the diagonal part 
and the off diagonal part as, 6M'^{(/))ij = M'^{(f))ij — M^{(/))ii6ij. The divergent part of Viioop 
can be easily computed by expanding it up to the second order of 5M^, 



^Uoop = V^''^ + K, 

V^'^ = ^/^TrLn{(Dr + M|(0))5.,+5M^.-aim?,} 



«,j=i 
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where, 



Ml + ml^-k'^ (1<^<4), 
Ml + ml,-k^ (5<z<8). 



(A2) 



The diagonal parts of the propagators are given as, 



Da 



Ml+mf^-k^ (1<«<4), 

m-t4^ (5<^<8). 



(A3) 



In the modified minimal subtraction scheme, Feynman integration is carried out with help 
of the well known formulae of dimensional regularization, 

and. 



^ 2 J {2'KYi ^^ ' 



GAn'^e 6Att'^ 



log- "2 



(A4) 



/^ 



A-d 



cf^k 



(27r)'^i(m2-P)(m2-fc2) 



1 1 



div. 



167r2e' 



with I = i — log47r and e = 2 — ^. The divergent part of y*^^^ is. 



8 



i=5 



-^ ^ ((5M2 - ml^a{),j{5M^ - m\^a^),,, 



QAn^e 



i^j=i 



327r2e 



m?,^M2(<^) + m^2^M^(0) + 2K,+ 



m, 



22^ 



i=l 



j=5 



- ^Tr[(M2(0)-mVi)(M2(0)-"^?2^i) 

^Tr[M^]. 



(A5) 



(A6) 



The trace of Eq. (1A6l) is calculated in Eq. (IB6l) and Eg. (IB 11 II of appendix B and the result 
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IS, 



V, 



(1) 



1 



div. 



327r2e 



m2,{6Ai($I$i) + 2(2A3 + \a){^\^2)} 



+ ml^{2{2\^ + A4)($l$i) + 6A2($|$2)} 



+ 



2m?2 



(2A3 + 4A4)($I$2 + $^$i: 



647r2e 

8m^2 + 4(mf;^ + m\2 



GA-K^e 



{12X1 + 4A3A4 + 4A^ + 2A^)(<l>l$i^ 



2Vrf>t, 



+ (12A^ + 4A3A4 + 4A^ + 2A^)(<I>T<|.2)' 

+ (I2A1A3 + 4A1A4 + 8A^ + 4A^ + I2A2A3 + 4A2A4)($l$i)($^$2 



+ (4A1A4 + I6A3A4 + 8X1 + 4A2A4) I $I$2 I' 



(A7) 



Now the counterterms for the one loop effective potential are simply given by changing the 
sign of the divergent part of Eq. (IA7I) , 



Vr, = -V,. 



(1) 



div. 



{Ii[M^ 



(A8) 



Using Eq. OASP and Eq. (lA4|) . one can derive the finite part of the 1 loop effective potential 
given in Eq.fl?! 



Appendix B: Derivation of Eq. (lA7|) 

In this section, we present the derivation of Eq.( lA7|) . We start with the quartic interaction 
terms of the Higgs potential, 



y(4) 



Ai 



Aa 



E' 



i=l 



.¥ + ^(E*')' + t(E^.')(E' 



j=5 



i=l 



i=5 



+ -^ ((01<^5 + 02</>6 + (/"S^T + </'408)^ + {Mq + </'3</'8 " </'2</'5 " (t>'ihf) ■ (Bl) 
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By taking the derivatives of V^^\ one can obtain the mass squared matrix M^ ((/)). One first 
computes the first derivative of V^^^ with respect to (pi, 



d(j)i 



f 2(E;=i '^,')20. + f<P^ E,=5 ^] + WM5 + Me + 0307 + 0408)0.+4 
+ (0106 + 0308 - 0205 " 0407)('^H06 " ^2*05 + 53i08 " (^4107)}, (1 < ^ < 4) 



A2 ' 



A4, 



2(e;=5 0,)20. + f 0. e;=i r, + 1 {(0 



1^5 -r (P2(P6 



+ 



^3</^7 



+ 



^i-4 



+ (0106 + 0308 - 0205 " 0407)(-55i02 + 56J01 " <^7i04 + <^8i03)}- (5 < i < 8). 



(B2) 



The second derivatives are given as, 



^i^'Vj 



\ {5^^ ELl 0^ + 20,0.) + f ^.,(EL5 0D + T {0.+40.+4+ 

(5lj06 - '^2j05 + '^3j08 - 54j07)(^H06 - (^2105 + hi4>& " hih)} , {^ < h j < 4), 
A30i0j + ^{0j+40i-4 + Ylt=l 5i+4j0fe0fc+4 + {S^jh + kj4>l " (^7j04 + ^8^03) 
(5h06 - (^2i05 + (^3*08 " <^4i07) + (0106 + 0308 - 0205 " 0407) 

{Su6(ij + 63i6sj - S2i6r,j - 64iS7j)}, (1 < i < 4, 5 < j < 8), 

A30i0i + ^{0i-40i+4 + Z]fc=l (^i-4j0fe0fc+4 + (5li06 " '^2i05 + '^3^08 " 54i07)x 
(-<^5i02 + <^6i01 - <^7i04 + <^8i03) + (0106 + 0308 " 0205 " 0407) X 

{6u6qj + 63i6sj - S2iSr,j - 64i67j)}, {5 < i <8,1 < j < 4), 

f {S., EU 0^ + 20,0.) + ^UEt=l 4>l) + ^{0,-40.-4 + 

(-(^5i02 + (^6i01 - (57,04 + (58i03)(-'55i02 + (^6.01 " <57i04 + <58i03)}, (5 < i,i < ( 



With Eq.f lBSp . the diagonal sums of M^ are given as, 



(B3) 



^M2 = 3Ai ^02 ^ 2A3 J]0,' + A4^0? = 6Ai<l>l$i + (4A3 + 2A4)$^$2, (1 < ^ < 4), 



i=l 



i=l 
8 



i=5 
4 



i=5 

4 



^M^ = 3A2^0,' + 2A3X;0' + ^4 5^0? = 6A2$|<l'2 + (4A3 + 2A4)$I$i,(5<2<8). 

i=5 i=5 j=l «=1 

(B4) 
The counterterm in Eq. (IA8p includes the following contribution, 

Tr[(M2(0) - m?2^i)(M2(0) - ml^ai)] = Ti[M\(j)) M\(P) - 2ml2aiM^] + 8ml^. (B5) 
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The second term of Eq. (IBSP is proportional to, 

TrKsfTiM^] = {2X3 + AX4){M5 + Me + h<f>7 + Ms)ml2 

= (2A3 + 4A4)($I$2 + $^$i)m22. (B6) 

The first term of Eq. fIBSp can be decomposed as, 

4 4 8 

i,j=l «=1 i=5 

8 

+ ^ M2(0),,m2(0),,. (B7) 

Each term of Eq. flB7p is given as. 



2 



4 



5^ M^(0).,M^(0),, = 3XliY: ct>] + 3A,A3 5: <P^ Yl ^? + ^1^^ E '^^ E 



/,2 



+ (</'l05 + 0206 + 0307 + 0408)^ + (0106 + 0308 " 0205 " 040?)^} 

E^? +-^MS^? +1 E' 

j=5 / \j=5 / \j=5 

= 12A2($1$i)' + (I2A1A3 + 4AiA4)($I$i)($5$2) 

+4A1A4 I $I$2 r +(4A3A4 + 4A^ + 2A^)($^<I>2)', (B8) 

4 8 8 4 I' 4 4 

^ Y, M\ct>\,M\<t>),. = >^IY. <^' E '^l + 2-^3A4 <^ 5^ 0.0i+4 Y. ^3^+^ 

1=1 j=5 i=5 j=l (_ i=l j=l 

+ (0106 - 0205 + 0308 " 040?)^} 

^2 f 4 8 / 4 \ 2 

+ Y 1 E '^^^ E '^? + 2 E <^*'^*+4 + 2(0106 - 0205 + 0308 - 0407)' 

[ i=l j=5 V i=l / 

= (4A^ + 2X1) ($I$i)($^$2) + (8A3A4 + 4A^) I $1$2 r, (B9) 

5^ M^(0).,M^(0),. = 3A^ 5^ 0,^ + 3A2A3 Y ^^ E '^J + ^2A4 E '^^ E '^^ 

i,j=5 \i=5 / i=5 j=l L j=l j=5 

+ 0206 + 0307 + 0408)' + (0106 + 0308 " 0205 " 0407)'} 



''l'A'5 



+A3AJE*?) +AifE^?') +lfE«' 



2 , 

, i=l / \ i=l / \ i=l 



= 12A2($^$2)' + (I2A2A3 + 4A2A4)($l$i)($^$2) 

+4A2A4 I $I$2 r +(4A3A4 + 4A^ + 2A^)($l$i)l (BIO) 
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From Eq.dHS]), Eq.(lB9]),and Eq. dBTOJ) . one obtains 



Tt[M\<P)M\(P)] = (12Af + 4A3A4 + 4A^ + 2A^)(<l>{<l> 



+ (12A2 + 4A3A4 + 4A2 + 2A2)($t$2)2 

+ (I2A1A3 + 4A1A4 + 8A^ + AXl + I2A2A3 + 4A2A4)($I$i)($^$2) 

+ (4A1A4 + I6A3A4 + 8A^ + 4A2A4) I $I$2 r . (Bll) 



Using Eq-dBH), Eq.(lB5]), Eq.(lB6]), and Eq-dHH]), one can derive Eq.(jMl). 



T 9Af 2 . 



Appendix C: [O^^Ojjj and L/ 



/ 



In this appendix, we show [0^^^0]jj and Ljj which are needed to calculate one loop 
corrections to the order parameters if} in Eq. (l29|) . [0'^^^0]jj (/ = 1, 2, 3, 4) are given as, 



dipi 

^) in W.r, rt9Qll \nT 

dipi 



/ H^i + ^2 + 6A3 - 2A4 - cos(4/3) (Ai + A2 - 2 (A3 + A4))) \ 

i (Ai + A2 + 6A3 - 2A4 - cos(4/3) (Ai + A2 - 2 (A3 + A4))) 

I (Ai + A2 + 6A3 + 6A4 - cos(4/3) (Ai + A2 - 2 (A3 + A4))) 

12{A2 cos^ 7 sin^ /3 + cos^ /3 sin^ 7A1} + (3 cos 2(/3 - 7) - cos 2(/3 + 7) + 2) (A3 + A4) 

\l2{Ai cos^ f3 cos^ 7 + sin^ /3 sin^ 7A2} + (-3 cos 2(/3 - 7) + cos 2(/3 + 7) + 2) (A3 + A4)/ 

(C2) 



and. 



L a/3 J^^ 2 

A2 cos2(/3) - sin2(/3)Ai - cos(2/3) (A3 + A4) \ 

A2 cos\l3) - sin2(/3)Ai - cos(2/3) (A3 + A4) 

A2 cos2(^) - sin2(/3)Ai - cos(2/3) (A3 + A4) 

3A2 cos2(7) - 3 sin^ 7A1 + 2^(sin(2(/3 + 7)) - 3 sin(2(/3 - 7))) (A3 + A4) 

V-3Ai cos2(7) + 3 sin2(7)A2 - 2-i2^(sin(2(/3 + 7)) - 3 sin(2(/3 - 7))) (A3 + A4)/ 

(C3) 
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Next we show Ljj in Eq. (l30|) . Note that Lu is symmetric Lu = Ljj and its non-zero 
elements are, 



I/ii = cos^ /^mli + sin^ /3m22 — 2 cos(/3) sin(/3)m 



12 



1 

+ 2 



L 



22 



^12 



3i;2{Ai cos^(/3) + sin2(/3) (2 (A3 + A4) cos2(/3) + sm^i(3)\2) } 

(A2 - Xi)v 



9 , cos 4/3 , , ^ , , . X X 9 cos 2/3 , , , . o 

^{ 7^(Ai + A2-2(A3 + A4))t;' + — ^^^ ' '-^ 



+ 2m^2 sin 2/3 — cos 2/3(m^]^ — 777-22)}, 

sin 4/3 1 

= L21 = i;{ ^ (Ai + A2 - 2 (As + A4)) t;' + - sin 2/3 (A2 - Ai) v\ 

— 2m^2 cos 2/3 — sin 2/3 (m^]^ — 777-22)} 
1 



i^33 = --^'sin(2/3)(t>2sin(2/3)A4-4m22) 
L44 = ti^cos(/3) sin(/3)mi 



^12- 



(C4) 



Appendix D: Orthogonal matrix O in Eq. (|3ip 



we show the orthogonal matrix 


in Eq.(|3iJ). 
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— sin/3 














cos/3 
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— sin/3 
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sin 7 
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— sin/3 

















sin/3 
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— sin 7 













^ 





cos/3 














sin/3y 
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